Field theory simulation of Abelian-Higgs cosmic string cusps 
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We have performed a lattice field theory simulation of cusps in Abelian-Higgs cosmic strings. The 
results are in accord with the theory that the portion of the strings which overlaps near the cusp is 
released as radiation. The radius of the string cores which must touch to produce the evaporation 
is approximately r = 1 in natural units. In general, the modifications to the string shape due to the 
cusp may produce many cusps later in the evolution of a string loop, but these later cusps will be 
much smaller in magnitude and more closely resemble kinks. 
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I. INTRODUCTION 

Cosmic strings arise naturally in spontaneous symme- 
try breaking phase transitions in which the vacuum man- 
ifold after the transition is not simply connected Q . (For 
reviews see (2]||). A simple field theory which possesses 
such defects is the Abelian-Higgs model, 



1 



-v 2 ) 2 



with 



D^<p = {dp - ieA^cj). 
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We can choose units such that 77 = 1 and e = 1. We 
will also work in the "critical coupling" regime in which 
j3 = A/ (2e 2 ) = 1 so that in our units A = 2. 

The fields of an infinitely long, static, straight string 
are given by B 



J(r) 

A H = - 



e lti f{r) 
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= a{0) — and f{r) — ► 1 and a{r) — > 1 
The exact form of / and a must be found 



where /(0) = 
as r — ► 00. 
numerically. 

If the string core radius is small as compared to the 
radius of curvature of a string, the action can be approx- 
imated by the Nambu action, 



S = -V-f d 2 (V- det 



7 



(5) 



where [i is the energy per unit string length, 7 is the 
metric on the world sheet of the string and C, 1 , C 2 denote 
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the world sheet coordinates. We can choose parameters 
a and t in such a way that the position of the string is 
given at any time t by a function x(cr, t) that satisfies 



|x'(M)| 2 + |*(M)| 5 
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x'(cr,i) •x(er,t) = 0. 



The equation of motion is 



'{a,t)=x{a,t), 



(6a) 
(6b) 
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where x' denotes differentiation with respect to a and x 
denotes differentiation with respect to t. Equation ( |6b| ) 
means that the velocity of the string is always in the plane 
perpendicular to the string while Eq. (|6a| ) parameterizes 
the string with constant energy per unit a. The general 
solution of Eqs. @ and (]?]) is 



xfot) = -(a(cr 



t)+b{a + t)), 



(8) 



where a and b are arbitrary functions that satisfy |a'| = 
|b'| = 1. The functions a and b describe waves traveling 
in the positive and negative a directions, respectively. 

As long as the string remains smooth, its entire evolu- 
tion is given by Eq. (g) and the state of the string at all 
times is fixed by the functions a' and b'. The values of 
these functions are unit vectors, so we can consider them 
to trace out paths on the surface of the unit sphere. For 
a string loop, all quantities will be periodic in a, so the 
paths on the unit sphere will be closed. Furthermore, 
since a and b are themselves periodic, 



J daa!{a) = J dab' {a) = 0, 



(9) 



which is to say that the center of gravity of the paths of a' 
and b' is at the center of the sphere. Thus, for example, 
a' or b' could not lie entirely in a single hemisphere, 
and, in general, we would expect that there would be 
one or more crossings between the paths of a' and b'. 
By the same argument, there will be places where a' and 
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— b' have the same value, that is to say values of a a 
and Ob at which a'(a a ) + h'(ab) = 0. As a result, at 
time t = (crb — u a )/2 and position a — (pb + cr a )/2 we 
will find x'(a,t) = (b'(er + t) + a! {a - t))/2 = and 
\x(a,t)\ = \(b'(a + t) -a'(a -t)\/2= 1. Such a point is 
a cusp. The string there is momentarily moving (in the 
Nambu-Goto approximation) at the speed of light. 

As the string evolves toward the time of the cusp, the 
Nambu-Goto approximation will break down. Since x' = 
at the cusp, the segments of string point in the same 
direction leaving the cusp, and so there is the potential 
for overlap between the segments on the two sides. When 
the string cores overlap, the topological constraints which 
stabilize the string no longer operate, and the energy in 
the string can be released as radiation. This radiation is 
of cosmological interest because it could potentially lead 
to observable cosmic rays, and because it affects the rate 
at which oscillating loops lose their energy [5|~[lC|]. 

At some point when the cusp-related evolution is com- 
plete, the string will once again be smooth and non- 
overlapping, and will resume Nambu-Goto evolution with 
new functions a' and b'. One can consider the evolution 
at the time of the cusp to be a scattering process which 
takes an initial state of a' and b' and produces a final 
state with new a' and b' and a certain amount of outgo- 
ing radiation. 

To determine the amount of radiation emitted, and the 
resulting configuration of the string, we have written a 
lattice field theory simulation of the Abelian-Higgs field 
theory. 



II. THEORETICAL EXPECTATIONS 

The usual expectation is that the string can participate 
in cusp evaporation when the "cores" of the strings on 
the two sides of the cusp are overlapping. A parameter 
r gives the "core radius," i.e., the point at which overlap 
dynamics can begin, as shown in Fig. [|. We expect this 
radius to be on the order of = 1 in our units. 

In @ we showed that a generic cusp can be ana- 
lyzed in a frame in which Xg" is parallel to xo, where 
the subscript denotes quantities evaluated at the cusp. 
Throughout this paper we will work in such a reference 
frame. 

In section III of we calculated the amount of over- 
lap between the two branches of a string at a cusp. The 
result was that the two segments of the string overlap 
from the cusp out to a distance a c given by 




(10) 



We expect an amount of string of total length 2a c to be 
replaced by a short bridging string of length 2r or less, 
so that the energy released in radiation will be approxi- 
mately 2\x<j c . 




FIG. 1. Schematic shape for a cusp in a cosmic string. The 
shaded region corresponds to the core of the string. To the 
left of the vertical line, the cores overlap, and the energy could 
be released as radiation. The string appears thinner near the 
cusp because of Lorentz contraction. 



In a string of cosmological size, a c would be very much 
larger than r, and the approximation in Eq. ( [To|) would 
be very accurate. However, we cannot simulate a string 
with such a large ratio of length scales. (Our largest sim- 
ulation has tJ c /r ~ 17.) To improve the accuracy of our 
prediction, we proceed as follows. We assume that the 
string has Nambu-Goto evolution until the time of the 
cusp (which is confirmed by simulation, see below), and 
compute the positions <r+ and er_ on the two sides of the 
cusp at which the (Lorentz contracted) string cores just 
touch. It is possible that a + and <r_ are not uniquely de- 
fined by this condition, since perhaps a + could be moved 
away from the cusp while cr_ was moved toward the cusp 
to produce different points of overlap. In this case, we 
take the points which lead to the largest emission of en- 
ergy. We then assume that the string segment from er_ to 
cr + will be replaced by a bridging string at rest between 
x(er_) and x(cr+), to arrive at an estimate of the energy 
release. Since we do not know the radius r at which 
overlap dynamics become possible, we can consider the 
simulation to provide us with a value for this radius. 

How should the paths of a' and b' be affected by the 
evaporation of the cusp? Near the cusp, a'(er) ~ a = 
— xo and b'(cr) ~ h' — xo. Outside the region of overlap, 
they will retain these values, but inside that region we 
will have a bridging string. The bridging string will point 
in the direction of Xg" because that is the direction in 
which the two arms of the cusp spread apart. Since the 
motion of the string near the cusp is approximately Xo, 
which is parallel to x " , the bridging string will be created 
approximately stationary. Thus on the bridge, a'(er) = 
b'(cr) oc x " with a positive constant of proportionality. 
There are thus two cases, depending on whether x " is 
parallel or antiparallel to xo. If they are parallel, then 
b'(cr) on the bridge will have roughly the same value that 
it has on nearby unaffected portions of the string, but 
a'(a) on the bridge will have nearly the opposite value as 



2 



FIG. 2. Schematic of the evolution of the functions a' (solid 
arrows) and b' (dashed arrows) before (dashed line) and after 
(solid line) the cusp evaporation process. 

it does nearby, as shown in Fig. 0. Thus in this case we 
expect the path of a'(er) to be broken near the point of 
crossing and reconnected around the opposite side of the 
unit sphere. In the case that Xq" and x are antiparallcl, 
we will have the reverse behavior. The path of a' (a) will 
be pretty much unaffected, while b'(er) will be broken 
and reconnected. This effect of the cusp on the evolution 
of the functions a' and b' on the unit sphere is shown in 
Figs. | and § 

III. SIMULATION 

A. Lattice Action 

To put Eq. (0) on a lattice, we have used a somewhat 
different approach from those that have been used in the 
past [Tl|-|i"3|] . We put the system on a lattice in both space 
and time, with the scalar field 4> stored on the vertices of 
the lattice, the vector field on the links, and the field 
strength F^ v on the faces. The action can be specified as 
follows. 

Let s denote a lattice site, / a directed link from site 
Z_ to site l + , and / a directed face made of links f\, /2, 
/3, and fi. Let J2 S denote summation over all sites, 
summation over all links, counting each link once, and 
summation over all faces, counting each face once. 
Let Ylm — s denote summation over all links that end in 
site s and X)/|/ 1= ; denote summation over all faces that 
border link I. Then the lattice action is 




FIG. 3. The paths of a' and b' before the cusp. The 
dashed segments are the parts that will be affected by the 
cusp. The dotted curves are the paths around the back of the 
sphere. 




FIG. 4. The paths of a' and — b' after the cusp. One of 
the paths has been broken and reconnected around the back 
of the sphere. 



where 



C s = -V(cf> s )=-^X{\^\ 2 - V 2 ) 2 

Li = giDi<j>Di4> 
1 

T 



£f = —9fFf 
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Here Ax and At are the spatial and temporal lattice 
spacing; gi is the metric, gi = +1 if I is timclike, and 
gi = — 1 if I is spacelike; and g/ — gf x gf 2 , i.e., <?/ = +1 
if all sides of / are spacelike, and gf = — 1 if / has a 
timelike side. The discrete covariant derivative is given 
by 



D l <t> = A l { 
D l 4> = A l { 



e iA ' eA ><f>i_) 
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where A; denotes Ax if I is spacelike or At if / is timelike. 
The field strength is given by 

F f = A h A h ( A h A fi + A h A h + A/3% + A /4 %) . 

(16) 

The action of Eq. (^) is invariant under a discrete 
gauge transformation given by a number A s at each lat- 
tice site, 



A i + A i - 1 (A i+ -A i J. 



(17a) 
(17b) 



We will work in temporal gauge, Ai = for all timelike 
links /, which greatly simplifies the calculation. 

The equations of motion resulting from Eq. ( pd|) are 



E giA i 







dV 
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E A h9fF}=ji 
f\h=l 



(18a) 
(18b) 



where the link current is 



jl = ieA [ 1 g t (e 



-ieA,A, 



eA,A 



(19) 



The stability of our finite difference scheme requires 
the Courant condition, At/ Ax < 1/^3 « 0.58 fj, for 
three spatial dimensions. In fact, ratios slightly less than 
this are still unstable in our code, presumably due to 
numerical error. To be safe, we have used At/ Ax = 0.5. 
To verify that this condition is sufficient, we have checked 
that our simulations conserve energy over long periods of 
time. 

We have used reflecting boundary conditions (fixed 
field values) on the top and the bottom faces of our lat- 
tice box, and first order absorbing boundary conditions 
| p5| , p^ | on the lateral faces. The absorbing conditions 
prevent the radiation emitted at the cusp from reflecting 
off the walls and interacting with the later evolution of 
the string. 



B. Initial Conditions 

To simulate cusp formation we must have an initial 
field configuration which will evolve into a cusp. This can 
be done exactly, using a result of Vachaspati et al. |l7]] 
that a traveling wave along a straight string is an exact 
solution of the field theory equations of motion. In other 
words, an arbitrarily shaped wiggle that propagates at 
the speed of light along a straight string will not radiate 
and will retain its form. 

These wiggles can be seen also in terms of functions 
a' and b' in the Nambu-Goto approximation. For exam- 
ple, a straight part of the string in the z direction can 
be written in terms of a' = z and b' = z. Introducing 
a wiggle traveling in one direction along the string will 
modify one of the functions a' or b' depending on the di- 
rection of the movement of the wiggle, keeping the other 
function fixed. 

Using this method we can obtain the Nambu-Goto 
description for two wiggles traveling toward each other 
along a straight string. Then by using [ fl7[ we can write 
down the expression for the fields in terms of the func- 
tions a' and b' and the static string fields given by Eqs. 
(0) and . By calculating the values for the fields when 
the two wiggles are still separated by a portion of straight 
string we ensure that our initial conditions are exact, so 
we do not have to use a relaxation procedure as is neces- 
sary in other field theory simulation schemes |Tl|-|T3]] . 

We design the two colliding wiggles such that they will 
produce a cusp in the Nambu-Goto evolution and use 
them as our initial conditions in the field theory code. 
One could imagine that the string would evolve into a 
different configuration due to interaction before the time 
of the cusp, but as we will show later, that does not oc- 
cur: the string core in the field theory simulation follows 
the path that the Nambu-Goto evolution predicts until 
immediately before the cusp. 

As described earlier, the cusp is characterized in the 
Nambu-Goto approximation by the intersection of a' and 
— b' on the unit sphere. Far outside the area of interac- 
tion, the string will be straight with a' at the north pole 
of the unit sphere, and — b' at the south pole. Thus we 
are looking for functions a' and b' which make loops on 
the sphere, starting and finishing at the poles and cross- 
ing at some point of the sphere. 

Two closed loops cannot intersect in just a single point, 
unless it is a point of tangency, which does not give rise 
to a generic cusp. But we want to obtain information 
about the energy radiated by just one cusp, so we do 
not want to have any interference in the dynamics due 
to the formation of a second one. (A realistic cosmologi- 
cal string would have multiple cusps, but they would be 
very well separated.) To prevent this interference wc will 
design the two wiggles so that the first cusp will alter the 
string in such a way that the second one will not occur. 
We will choose our parameters so that Xq" • xo > 0, so we 
expect a' to be broken and reconnected, while the path 
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of b' will not be greatly affected. Thus we will arrange 
for b' to cross a' twice at the same point, so that when 
a' is reconnected to eliminate the first crossing, the sec- 
ond crossing will be eliminated as well. We will choose 
a simple path for b' that starts at the pole, runs along 
a meridian, crosses a' at the equator, and continues for 
short distance before turning around and returning along 
the same path. 

We must also arrange that the magnitudes of a ' and bg 
are different to avoid a pathological situation which leads 
to extra energy emission, as discussed in fnl. We also 
want to work in the canonical frame where Xg is parallel 
to Xo, so we require that a " + b " is perpendicular to a ' 
and to bg. 

We will start by specifying our desired values for a ' 
and a ". As much as possible we would like other deriva- 
tives not to contribute near the cusp, so we will set 



a unnormalizcd( cr ) ~ a + aa + ~2 &0 



and then 



^unnormalizcd V^) 



(20) 



(21) 



unnormalizC' 



We will use this form for a € [— a\, oi], where o\ is chosen 
to be large enough that changes in a' outside this range 
will not have too much effect on the cusp dynamics. 

Outside of this range, we just have to connect a' to 
the north pole of the unit sphere. We do this by choos- 
ing paths of circles on the sphere which will smoothly 
interpolate between a'(— a{) and a' '(<J\) and the pole. 

The result of this design for a' and — b' is shown in 
Fig. |. 

To produce the initial field configuration, we offset the 
wiggles produced by a' and b', so that one is at the top of 
the lattice moving downward and other is at the bottom 
of the lattice moving upward. The actual shape of the 
string can be seen in Fig. ^| 

When the two wiggles combine, we will look at the 
shape of the string at the time of the cusp in the Nambu- 
Goto approximation, and define the "length scale" of the 
cusp for the sake of discussion as 



L = 



(22) 



C. Recovering the energy and form of the string. 

In order to obtain the amount of energy released in the 
cusp evaporation, we implemented a method to recover, 
out of the field theory simulation, the total amount of 
energy left in the string. 

We first identify the plaquettes of the lattice that the 
string goes through, by looking for places where the phase 
of the complex scalar field on the four corners of the 




FIG. 5. Initial conditions for the functions a' and — b', on 
the unit sphere. 




FIG. 6. The initial conditions. The surface shown is a sur- 
face of constant energy density surrounding the string core. 
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plaquette wraps around the origin in the complex plane. 
As long as two strings do not go through the same lattice 
cube, this technique always allows the string to be traced 
from one face to the next. We take the path of the string 
to go through the centers of these plaquettes. (We have 
tried some algorithms for improving the path of the string 
by closer examination of the field values at the corners 
of the plaquettes, but they did not produce significantly 
better results.) 

This technique produces a discrete approximation to 
the string path, which we then smooth by a Gaussian 
convolution. The end result is a list of smoothed strings 
at different times, parameterized by an arbitrary param- 
eter A. We will denote these string positions by x(X,t). 

Near the time and position of the cusp, the dynam- 
ics are not well described by the motion of a string, but 
rather involve the complete field configuration. One can 
still recover a path of faces around which the field phase 
winds, but the different positions of this path at different 
times do not represent motion, but rather the unraveling 
of a field configuration. However, since we are only inter- 
ested in the difference between the energy in the string 
before and after the cusp, we will just apply this method 
at times well before or well after the cusp where we are 
certain that the Nambu approximation is a faithful de- 
scription of the string. 

To recover the energy of the string we first calculate 
x(A, t) by computing 



x(A,t) 



x(A, t) - x(A', t - NAt) 
NAt 



(23) 



where N is a number of lattice time steps over which 
we average to reduce the effect of errors in recovering 
the string position, and A' is determined by the condi- 
tion that (x(A, f) — x(A',t— NAt)) be perpendicular to 
dx(X,i)/dX. 

Once we have calculated the velocity we can obtain the 
relation between A and the usual parameter a by using 
Eq. (§3), 



dx(A, t) 



dX 



da = 



dX 
da 

dX 



Vl-|x(A,t)|^ 



dx(X,t) 



Vi-|*(V)| s 



dX 



(24) 



(25) 



The total length of string is just the integral of the 
equation above and so the energy is 



Estring = \l \ d(J . (26) 



Once we have obtained the vectors x' and x for the 
string, we can get the functions, a'(<r, t), and h'(a, t) from 
the relations 



a'(a,t) = x'(a,t) - x(a,t) 
b'(a,t) = x.'(a,t) + ±(a,t) 



(27a) 
(27b) 



and study the change of these functions due to the cusp 
dynamics. 



IV. RESULTS 

We ran the simulation code for cusp scales (as defined 
by Eq. (|22|)) from L w 5 up to I « 34. For sizes up 
to L sa 17 we used Ax — 0.25 or less, which gives a 
quite accurate simulation. For larger sizes we were forced 
by memory constraints to use larger values of Ax up to 
Ax = 0.425. These larger values produce somewhat de- 
creased accuracy, especially in those parts of the string 
which are rapidly moving. They also reduce the accuracy 
with which the position of the string and thus the energy 
emission can be recovered from the simulation. However, 
we feel that these simulations are still sufficiently accu- 
rate to draw some general conclusions about the cusp 
evaporation process. 

The code was written in Lisp and executed on a Digital 
Alpha processor. The largest lattices required about 700 
MB of memory and one week of CPU time. 

We first run the code letting only one wiggle evolve 
on the straight string, setting b' = z. This experiment 
confirms Vachaspati's result |l7| , showing the traveling 
wave on the string without any radiation. This is also 
a check for the code, since it shows stability and energy 
conservation for the evolution with the boundary condi- 
tions specified above. 

With the two colliding wiggles, the string forms a cusp 
as expected. Energy is conserved to within 1% until the 
radiation emitted at the cusp begins to be absorbed at 
the boundary. Figures |^ and || show the string shortly af- 
ter the moment at which the Nambu-Goto approximation 
predicts the cusp. Except near the cusp, the Nambu- 
Goto evolution has been followed accurately. But near 
the cusp, the core of the string has collapsed away from 
the area where the cusp energy is stored and has released 
this energy in the form of radiation. This is what one 
would expect from the overlap model. At some point 
near the cusp the strings are close enough that the fields 
can untangle, and the cusp beyond this point is replaced 
by a bridging string. After this time, the released radia- 
tion travels freely in an expanding shell as shown in Fig. 

For each simulation run, we recover the shape of the 
string, the form of a' and b', and the total length of 
string at each time step. Figure [l^ shows a typical graph 
of the recovered points. There is significant variation in 
the individual points, both random variation due to the 
way in which the string core cuts through the lattice, 
and systematic variation due to inaccuracies in recover- 
ing the energy that depend on the string shape. However, 
there is a clear difference between the string length be- 
fore and after the time of the cusp. (At times very near 
the cusp the energy cannot be recovered so there are no 
corresponding points on the plot.) 
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FIG. 7. String evolution at the cusp for L » 17. The core 
of the string in the field theory simulation is shown light, and 
the Nambu-Goto prediction dark. The "fog" shows areas of 
high energy density. The string core has collapsed away from 
the point of the cusp, leaving the energy behind in radiation. 




FIG. 8. Close-up view of Fig. @ 




FIG. 9. Close-up view of the core of the string and the 
radiation after the cusp evaporation. The cloud of radiated 
energy travels in an expanding shell. 
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FIG. 10. Recovery of the total amount of a on the string 
during a cusp evaporation, for L w 17 . Points at left are be- 
fore the cusp; points at right are those after cusp evaporation. 
Near the time of the cusp, a cannot be recovered. 
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FIG. 11. Amount of a as a function of the length scale 
L. The dashed line corresponds to the theoretical prediction 
assuming r = 1.15. 
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Radius 

FIG. 12. Functions / (solid curve) and a (dashed curve) 
for the static straight string. The vertical line indicates the 
best-fit value of r from the field theory simulation. 

In Fig. [ll] we plot the string length released by the 
cusp and compare with theoretical predictions. The er- 
ror bars in Fig. 
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are a combination of the variation in 
the actual recovered energy with an estimate of the sys- 
tematic errors due to the recovery process. As shown, the 
results are in accord with the overlap model with over- 
lap radius r ss 1.15. Figure |lj shows this value of r as 
compared to the shape of the string profile. One would 
expect the overlap dynamics to begin at a radius where 
the fields have values intermediate between the core val- 
ues and those at large distances, and the result seems 
quite reasonable in this regard. 

The recovered a' and — b' shortly after the cusp are 
shown in Fig. 13. By comparison with Fig. ||, we can see 
that the path of — b' is roughly the same before and after 
the cusp, but the path of a' has been greatly altered. Far 




FIG. 13. Recovered functions a' and — b' after the cusp, 
for L » 17. 



from the cusp it is the same, but in the region of inter- 
action, it has been reconnected around the other side of 
the sphere from the original crossing point. This recon- 
nection, as planned, prevents the occurrence of a second 
cusp. 



V. DISCUSSION 

We have simulated a cosmic string cusp in Abelian- 
Higgs field theory on the lattice. Our results are in accord 
with the model that the evolution is accurately given by 
the Nambu-Goto equations of motion until nearly the 
time of the cusp. The effect of the cusp is to release 
an amount of energy which is well approximated by the 
overlap model with r ~ 1.15. 

As discussed in for a string of cosmological size, 
this model gives an energy emission of order fi^/rL, where 
L is the typical length scale. This is smaller than previous 
estimates by a factor of (r/L) 1 / 6 and would not lead 
to observable cosmic rays. 

The effect of the evolution through the cusp is to dis- 
connect the path of either a' and — b' and reconnect it 
around the other side of the unit sphere. We will suppose 
without loss of generality that x " • xo > so that it is 
a' that is reconnected, while — b' merely has a section 
deleted and forms a small kink. It is likely that the re- 
connected a' will have many crossings with — b' and thus 
might lead to many cusps in the future. However, the 
total amount of string length in which a' loops around 
the sphere is quite small, on the order of r. Thus the re- 
maining feature involves a bending of the string in about 
a string radius, and is thus more like a kink than a re- 
gion of string which would produce a cusp. Any cusps 
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in which it is involved will have even smaller emission of 
radiation than the original cusp. 

This result could perhaps be modified by gravitational 
radiation, which smoothes out small-scale features p8|. 
Gravitational back-reaction might smooth out the region 
where a' varies rapidly and restore it to its original cross- 
ing with — b'. This would produce another cusp at the 
same position in a subsequent oscillation. However, the 
region of overlap contains much less energy than before, 
because most of it has been emitted, so that a' and — b' 
cross this region in a small amount of a. Subsequent 
gravitational effects might increase this energy, but it 
will not grow as large as in the original cusp, and so will 
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